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Abstract 



We study the phase transition of the ±J Heisenberg model in three di- 
mensions. Using a dynamical simulation method that removes a drift of the 
system, the existence of the spin-glass (SG) phase at low temperatures is 
suggested. The transition temperature is estimated to be Tsg ~ 0.18 J from 
both equilibrium and off-equilibrium Monte-Carlo simulations. Our result 
contradicts the chirality mechanism of the phase transition reported recently 
by Kawamura which claims that it is not the spins but the chiralities of the 
spins that are ordered in Heisenberg SG systems. 



Typeset using REVTeX 



1 



It has been believed that the spin-glass (SG) phase is realized in three dimensions (3D) 
for the Ising model [|T],[| but not for the XY and Heisenberg models [|3|-fJ. Thus, the SG 
phases observed in experiments were suggested to be realized due to magnetic anisotropies 
However, it remains puzzling that no sign of Heisenberg-Ising crossover has been 
detected, which is expected if the observed finite temperature phase transition is due to a 
weak magnetic anisotropy. Recently, Kawamura and coworkers have proposed a chirality 
mechanism to solve the puzzle [PHT2|. They considered the chirality described by neighboring 



three spins. It has either positive or negative value, just like the Ising spin. Kawamura fll 
calculated the chirality autocorrelation function C x (t) and the spin autocorrelation function 
Cs{t) of 3D Heisenberg SG models with and without a weak random magnetic anisotropy D 
and found that, even for D = 0, C x (t) exhibits a pronounced aging effect reminiscent of the 
one observed in the mean- field model, while Cs{t) exhibits a similar aging effect only when 
D^O. He claimed that, in the 3D Heisenberg SG model, a chiral-glass (CG) phase transition 
occurs at a finite temperature T CG 7^ 0, but the SG phase is absent. He argued that, in real 
SG magnets, the spin and the chirality are mixed due to a weak magnetic anisotropy, and the 
CG transition is revealed as anomalies in experimentally accessible quantities. According 
to this interpretation, the SG phase transition never occurs in Heisenberg-like systems, and 
what was observed in experiments is nothing but the CG phase transition. This mechanism 
is quite interesting, because it calls for re- consideration of the SG phase transition from 
both theoretical and experimental points of view, because most of SG magnets are rather 
Heisenberg-like. 

In this Letter, we demonstrate that the apparent difference between Cs{t) and C x (t) 
comes from a drift of the system and suggest that, contrary to the chirality mechanism, the 
SG phase is realized in the Heisenberg SG model even when the anisotropy is absent. We 
estimate the SG transition temperature of the ±J Heisenberg model as T$q ~ 0.18J from 
both equilibrium and off-equilibrium Monte-Carlo simulations. We believe that what was 
revealed by Kawamura and coworkers is nothing but the ordering of the spins themselves. 
The most important point revealed here is that the SG phase will be much more stable than 



what were believed so far and re-examination is also necessary for different models |13|-[T8|. 

We start with the ±J Heisenberg model on a simple cubic lattice of L x L x (L + 1)(= N) 
described by 

H = -^J ij S i S j , (1) 

(ij) 

where Si is the Heisenberg spin of \Si\ — 1 and (ij) runs over all nearest-neighbor pairs. 
The exchange interaction (bond) Jjj takes on either + J or — J with the same probability of 
1/2. 

First we consider the dynamical properties of the system calculating the spin autocorre- 
lation function. When one considers the dynamics of the Heisenberg model with D = 0, one 
should pay a special attention to a drift of the system, i.e., the global rotation of the system. 
Although this drift does not affect the static spin correlation and becomes negligible in large 
systems, some data of finite systems inevitably suffer from this effect. The spin autocorre- 
lation Cs{t) diminishes by this drift, whereas the chirality autocorrelation C x (t) does not, 
because the chiralities are invariant under the global rotation of the system. Therefore the 
difference in the behavior between Cs{t) and C x (t) might come from this drift. 

In order to consider the dynamical properties of an infinitely large system, we remove 
this effect applying the uniform rotation to the system so that 

S = y £\R{t)S i (t + t v) ) - S^ 2 (2) 

i 

becomes minimum, where R(t) is the rotation matrix, and {Si(t w )} and {Si(t + t w )} are 
the spin configurations at times t w and t + t w , respectively. We can successively determine 
R(t), because R(t) = E at t = and it changes step by step as t goes by, where E is the 
unit matrix ||19|| . Hereafter, we call the system described by {Si(t)} the drifting system and 
that by {Si(t + t w ) = R(t)Si(t + t w )} the fixed system. The spin autocorrelation function 
Cs(t w , t + t w ) of the fixed system is defined as 

C s (t w , t + t w ) = — J2l(Si(t + t w )Si(t w ))], (3) 

i 
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where (• • •) and [• • •] mean the thermal average and the bond distribution average, respec- 
tively. The chirality autocorrelation function is also defined as 

C x {t w , t + t w ) = — Y^[(Xi»{t + t w ) X iM)], (4) 

where x%n is the chirality at the ith site and in the /ith direction defined by Xip = &i+e ' 
(Si x ) with (// = x,y,z) being a unit lattice vector along the \i axis. Note again 

that the chirality Xi^ is invariant under the global spin rotation, and then the chirality 
autocorrelation functions of both drifting and fixed systems are completely the same and 
described by Eq. (4). Before applying this method to the present model, we test it in the 
ferromagnetic Heisenberg model, results of which are shown in Fig. 1. It is found that 
the spin autocorrelation function of the fixed system with a large t w approaches to the 
equilibrium value of the square of the order parameter, ((M/N) 2 ), with M being the total 
magnetization, whereas that of the drifting system rapidly decays with increasing t. This 
result indicates that the effect of the drift can be removed properly by the present method. 

We perform the same simulation that was done by Kawamura JO]] using the standard 
single spin-flip heat-bath Monte-Carlo(MC) method. That is, starting with a random initial 
spin configuration, the system is quenched to a working temperature, and after waiting t w 
MC steps per spin (MCS) the autocorrelation functions are measured up to about 2 x 10 5 
MCS. A sample average is taken over about 32 independent bond distributions. The results 
of Cs(t w ,t + t w ) for different t w are presented in Fig. 2 as functions of t. Contrary to 
that reported by Kawamura, Cs{t w ,t + t w ) exhibits a quite similar aging effect that was 
found for C x (t w ,t + t w ) (see Fig. 1 in ref. fill), which indicates that the spins and the 
chiralities possess similar dynamical properties. To examine this point in more detail, we 
plot in Fig. 3 the ratio T(t w ,t + t w ) = C x /C$ for a fixed t w = 10 5 at different temperatures. 
At all the temperatures, T(t w ,t + t w ) slowly decreases as t goes by, revealing that Cs decays 
more slowly than C x . This result clearly indicates that, if the ordering of the chiralities is 
realized, then the same is true for the spin. We estimate values of the Edward-Anderson 
order parameters of the spin and the chirality, gsc and gcc, using the method proposed by 
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Parisi et al. |2(J and used by Kawamura |TT|. That is, these values are extracted by fitting 
the data of C$ and C x for t w = 3 x 10 5 to the power- law form of 

C 

Cs,x(tw,t + t w ) ~ gsG.CG + 7v (5) 

in the time ranges of 30 < t < 3000 and 50 < t < 5000 . The obtained gsc an d <7cg are 
plotted as functions of temperature in Fig. 4. Both quantities have positive, non-vanishing 
values at low temperatures and seem to vanish at almost the same temperature of T ~ 0.18 J 



21fl . This result suggests that, in fact, the phase transitions of the spin and the chirality 
exist at T ^ and they occur at the same temperature. It is not easy, however, to estimate 
the definite value of the transition temperature as well as values of the order parameter 



exponents |p2]j , because the data of gsc around T = 0.18 J depend strongly on the time 
range used for extracting them. 

Next we perform a detailed MC simulation in the thermal equilibrium to examine the 
occurrence of the SG phase transition itself. We calculate the SG susceptibility xsg defined 
by 

XsG = ^T,{(^S j n (6) 

ij 



using the exchange MC algorithm p3| . The sizes of the lattice studied here are L = 5 ~ 19. 
Equilibration is checked by monitoring the stability of the results against at least two-times 
longer runs. The numbers of the samples are 320 for L < 11, 96 for L = 15, and 48 for 
L = 19. The results of xsg f° r different L are plotted in Fig. 5. At low temperatures, xsg 
exhibits a strong size dependence implying the occurrence of the phase transition at T ^ 0. 
If the lower critical dimension di is less than the lattice dimension, di < 3, and the phase 
transition really occurs at T = T$q, the data for different L will be scaled as 

XSG = L 2 ^F(L^(T-T SG )), (7) 

where v is the exponent of the correlation length and r] is the exponent which describes 
the decay of the correlation function at T = T$q. The scaling plots obtained by assuming 



Tsg 7^ and T$g = are shown in Fig. 6. It should be evident the scaling with Tsg ^ 
works much better than that with Tsg = 0, even if the data for the smallest size L = 5 are 
ignored in the latter ||24| . The phase transition temperature and the values of the critical 
exponent estimated here are Tsg/ J = 0.18 ± 0.01, v = 0.97 ± 0.05 and r\ = —0.1 ± 0.1. 
Note that, in the ± J model with the random magnetic anisotropy D ^ 0, the transition 
temperature was estimated as Tsg ~ 0.30J for D = 0.1J, and Tsg ~ 0.25J for D = 0.05J 
||. Thus our value of Tsg ~ 0.18 J is consistent with the one obtained by extrapolation of 
those values. The value of v ~ 1.0 for the correlation length also is not strongly different 



from that of v ~ 1.2 estimated in the case of D ^ p5| . 

In summary, all the results presented here seem to support the view that the SG phase 
is realized in the ±J Heisenberg model and the phase transition to this phase occurs at 
Tsg ~ 0.18J. This view is compatible with our recent studies of the stiffness of the system, 
which suggest that, contrary to previous studies [^f|,|9],[l(J , the stiffness exponent 9s of the 
spin has the positive value of 9s ~ 0.8 at T = |^6[ and changes its sign at T ~ 0.19 J [|27| . 
Thus we believe that the SG phase is realized in this model even when the anisotropy is 
absent. A major role of the anisotropy will be neither to bring the phase transition nor to 
mix the spin and the chirality, but to fix the system and lead anomaly of the susceptibility in 
a finite system ||. One might think that the present suggestion is quite surprising, because 
it disproves the common belief. It should be noted again that the occurrence of the ordered 
phase of the chiralities was already suggested by Kawamura et al. and what is suggested here 
is that the spins themselves also order. We think the origin of the occurrence of the phase 
transition in the Heisenberg SG model is the ordering of the spins and that the ordering of 
the chiralities happens to be accompanied by it. However, further studies are necessary to 



establish this concept as well as to reveal the nature of the phase transition [28]. We hope 



that the present paper stimulates studies of SG systems with the continuous spin symmetry. 
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FIGURES 

FIG. 1. The spin autocorrelation functions of the ferromagnetic Heisenberg model in 3D 
plotted versus t for a waiting time t w = 2 x 10 4 at temperatures around the Curie temperature of 
Tq/J ~ 1.45. The filled and open symbols are data from the fixed and the drifting systems, respec- 
tively. The arrows indicate the equilibrium values of the square of the magnetization ((M/N) 2 }. 

FIG. 2. The spin autocorrelation functions of the ± J Heisenberg model plotted versus t for 
different waiting time t w . The data are averaged over 32 samples of the lattice size L = 15. 

FIG. 3. The ratio of the chirality to the spin autocorrelation functions T(t w ,t + t w ) = C x /Cs 
of the ± J Heisenberg model plotted versus t. The data are averaged over 32 samples of the lattice 
size L = 15. 

FIG. 4. Temperature dependences of the Edwards- Anderson order parameters of the spin and 
the chirality, qsc and qcc, of the ±J Heisenberg model. Open and solid symbols indicate those 
extracted in the time ranges of 30 < t < 3000 and 50 < t < 5000, respectively. The data are 
averaged over 32 - 200 samples. Lines with (3 = 0.5 for qsc and /3qg = 1-0 for qcc are guides to 
the eye. 

FIG. 5. Temperature dependences of the spin-glass susceptibility xsg of the ± J Heisenberg 
model in 3D for different sizes of the lattice. 

FIG. 6. Typical examples of the finite size scaling plot for (a) T$g ^ and (b) T$q = 0. 
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Fig. 1 F. Matsubata et al 



L = 15 T/J = 0.05 




4 i 1 1 1 1 1 1 1 1 1 1 1 1_ 

' 10 1 10 2 10 3 10 4 10 5 

t 



Fig. 2 F. Matsubara et al. 
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